Hamiltonian dynamics reveals the existence of quasi-stationary states for long-range 

systems in contact with a reservoir 
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We introduce a Hamiltonian dynamics for the description of long-range interacting systems in 
contact with a thermal bath (i.e., in the canonical ensemble). The dynamics confirms statistical me- 
chanics equilibrium predictions for the Hamiltonian Mean Field model and the equilibrium ensemble 
equivalence. We find that long-lasting quasi-stationary states persist in presence of the interaction 
with the environment. Our results indicate that quasi-stationary states are indeed reproducible in 
real physical experiments. 
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The statistical mechanics of systems with long-range 
interactions is important for a variety of physical appli- 
cations, including, e.g., gravitational systems, plasmas, 
Bose-Einstein condensates In such systems, the inter- 
particle interactions decay at large distances r as l/r" 
with a < d (spatial dimension) and ordinary statistical 
mechanics assumptions are questioned by nontrivial ef- 
fects, like persistence of correlations and non-negligible 
interface energies. In particular, the Boltzmann trans- 
port equation picture for the approach to equilibrium is 
not valid and long-range interacting systems may even 
display inequivalences among different equilibrium statis- 
tical ensembles 0, Q . Because of these subtleties a priv- 
ileged investigation tool is the microscopic Hamiltonian 
dynamical simulation. Hamiltonian dynamics at fixed- 
energy (microcanonical ensemble) for a paradigmatic 
long-range Hamiltonian (see below) directly connected 
with experiments revealed the existence of long-lived 
quasi-stationary states (QSS) that finally cross over to 
Boltzmann-Gibbs (BG) statistical equilibrium The 
question arises about the reproducibility of such QSSs in 
real physical experiments, where the environment intro- 
duces perturbations that cannot be taken into account 
in the microcanonical ensemble. If so, practical advan- 
tages could be obtained by knowing control mechanisms 
that improve or hinder this quasi-stationarity. The use 
of classical numericalprescriptions (like the Nose-Hoover 
or the Monte Carlo |3)i to perform such an investigation 
raises some subtle questions, since these methods implic- 
itly assume the BG equilibrium. Here we propose a novel, 
physically transparent, microscopic Hamiltonian dynam- 
ics obtained by coupling the long-range system with a 
thermal bath (TB) which simulates the effect of the envi- 
ronment on the system itself. For sufficiently large time- 
scales, such dynamics confirms the equilibrium canoni- 
cal ensemble predictions. However, starting with out-of- 
equilibrium initial conditions, we discover that long-lived 
QSSs subsist in presence of the TB, providing evidence 



in favor of the reproducibility of QSSs in real physical ex- 
periments. We also discuss the relaxation process follow- 
ing the QSS and the peculiar behavior of the Boltzmann' 
//-function. 

In a magnetic context, the Hamiltonian Mean Field 
(HMF) model H describes a set of M globally coupled 
Xy-spins with Hamiltonian 
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where Oi G [0, 27r) are the spin angles and li G M. their 
angular momenta (velocities). The presence of the ki- 
netic term naturally endows the system of spins with 
an Hamiltonian dynamics. This Hamiltonian is consid- 
ered "paradigmatic" for long-range interacting systems 
m since its equilibrium properties are analytically solv- 
able both in the microcanonical and in the canonical 
ensemble 0, Q and it is representative of the class of 
Hamiltonians on a one-dimensional lattice in which the 
potential is proportional to X]fj=i [1 ~ cos(6'i — 9j)] /rfj, 
where r,;,- is the lattice separation between spins and 
a < 1 (notice that the potential in Eq. ^ is re- 
covered in the limit a — * 0). Also, direct connections 
with the problem of disk galaxies 9] and free electron 
lasers experiments 5] have been established. Whereas it 
has been recently proven that Q does not present mi- 
crocanonical/canonical inequivalence at equilibrium 0], 
its unusual dynamical features received recently a lot of 
attention 0, M lldl fl^ . In fact, fixed-energy dynam- 
ical simulations starting with out-of-equilibrium initial 
conditions display the existence of longstanding (infinite- 
standing in the thermodynamic limit) QSSs appearing 
after a "violent relaxation" dynamics. During the QSS, 
phase functions such as the specific kinetic and potential 
energies fiuctuate around stationary or quasi-stationary 
non-equilibrium average values. In this Letter we intro- 
duce a microscopic setup where the HMF model is in con- 
tact with a short-range thermal bath (TB) in such a way 
that the thermodynamic limit is achieved with a negli- 
gible interaction energy. Equilibrium dynamics confirms 
the microcanonical/canonical equivalence for the HMF 



2 




FIG. 1: Sketch of the interactions considered in our canonical 
setup. Dashed lines mimic the interactions between the HMF- 
(fuU circles) and the TB- (empty circles) spins. Full (dotted) 
lines represent the HMF (TB) couplings. 

model. If the HMF-TB coupling is weak enough, the re- 
laxation to equilibrium is still characterized by drastic 
slowing-downs (QSSs) where also the system energy fluc- 
tuate around quasi-stationary average values. We discuss 
in a separate paper l^j the details about the statistical 
mechanics of QSSs in the canonical ensemble, a question 
of considerable debate 0, 0| . Here we report that 
in Gibbs' F-space their statistical mechanics is obtained 
using the classical BG definition of the entropy. 

The TB we consider is characterized by iV ^ Af equiv- 
alent spins first-neighbors coupled along a chain 

TV ,2 N 

Htb= i+ E [1 - cos(0,+i - 0,)] , (2) 

i=M+l i=M+l 

with 9n^i = 6m+i- The interaction between and (0) 
is modulated by a coupling constant e: 

M S 

i/, = e^^[l-cos(0,-0,^(,))] , (3) 

i=l s=l 

where rs{i) are independent integer random numbers in 
the interval [M + 1,N]. In this way, each HMF-spin is 
in contact with a set of S different TB-spins chosen ran- 
domly along the chain (see Fig. ^ . This set is specified as 
initial condition and remains fixed during the dynamics. 
The total Hamiltonian H = Hhmf + Htb + Hj defines 
then a microcanonical system (constant energy iJ), in 
which the energy of the HMF model can fluctuate. In 
our approach, the temperature is defined by (twice) the 
specific kinetic energy and we expect the TB to main- 
tain a constant temperature about which the HMF model 
thermally equilibrates. By assuming a "surface-like ef- 
fect" S ~ M'^-'^ (with < 7 < 1), we make sure that 
the interaction energy, Ej ~ Af*", satisfy Ehmf ^ M 
{Etb ^ N :s> M), thus ensuring a well defined thermo- 
dynamic limit. For the present results we chose = Af^ 
and S — 10^ M^^^^ . To integrate the equation of motion 
we use a velocity Verlet 7] algorithm with an integration 
step guaranteeing conservation of total energy within an 
uncertainty of AE/E ~ 10^^. 

Let chmf = Ehmf /M = [knMF + (1 ^ ^'hmf)I'A^ 
where kuMF = Y^fii ^If^^^ ^^d tuhmf = 




FIG. 2: (a): Time evolution of the Thmf and Ttb tem- 
peratures for A/ = 10^, e = 0.0 1 and Tp = 0.5. Initially, 
Phaif{1,0) = exp(-;V2r(t*))/v'27rr(t*) with T{t*) = 0.7. 
(b): Velocity PDF at t > t* . Solid line is ptb{1,0)- (c): 
Caloric curve. Solid line is the BG equilibrium and dashed 
line is the prolongation of the ordered phase to subcritical 
energies. Empty symbols are the average value of eHMF(t) 
at equilibrium. Stars refer to Nose-Hoover simulations. Full 
circles correspond to the QSS studied in the paper and to the 
microcanonical and canonical equilibrium obtained as t — > oo. 



\^i=i{cos9i,sin9i)\/M are respectively the specific 
kinetic energy and the magnetization of the system. 
It is known that at crmf = 0.75 and temperature 
Thmf = 0.5 (in natural dimensionless units) a contin- 
uous phase transition occurs separating a disordered 
{niHMF = 0) phase from a ferromagnetic one 0|. 
Here we show that our Hamiltonian dynamics con- 
firms such equilibrium predictions. The width Tq 
of the Maxwellian PDF for the initial TB-velocities 
Ptb{1,0) = exp{—P /2To) /^/2ttTo is a control parameter 
through which we set the TB temperature. In fact, after 
a transient relaxation {0 < t < t* ^ 100) the TB reaches 
its own equilibrium at the target temperature Tq (i.e. 
2kTB{t) ~ToWt>t*). Att = t*we then switch on the 
HMF-TB coupling, Hj, by setting e{t) =e* >OWt>t*. 
For e* = 0, Hi ~ 0, and the scheme reproduces the mi- 
crocanonical dynamics of the HMF. The setup was tested 
for many different initial conditions of the HMF model 
with 10^ < A/ < 10* and 0.005 < e* < 0.1. In aU cases, 
for t ^ t* , the system reaches the thermal equilibrium 
characterized by 2kHMF{t) — Tq (Fig[21i) a velocity PDF 
PHMFil,t) — Ptb{1,^) (Fig. 13)), and an equilibrium 
magnetization. The relaxation to equilibrium could last 
very long and typically occurs through a number of 
drastic slowing-downs during which the average value of 
Thmf is constant or almost constant (plateaux in Fig. 
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FIG. 3: Microcanonical (a) and canonical (b,c) QSSs in terms 
of 

iTi^HMF- The curves have been obtained by averaging over 
a number I of initial conditions following the probability dis- 
tributions described in the text. I varies from a maximum of 
20 for small M down to 5 for big M. (d): Log-log plot of 
as a function of the system size M. The QSS life-time has 
been defined as the time at which a phase function (Thmf for 
e* = and euMF for e* 7^ 0) changes 10% of its stationary 
average value. 



Similar effects were found in using a stochastic 
dynamics (see also j^] for a stochastic canonical version 
of the HMF, named Brownian mean field). By varying 
To, we obtain an estimate of the caloric curve in excellent 
agreement both with the BG equilibrium prediction and 
with the microcanonical simulation, even at the critical 
temperature and independently of e* (Fig. These 
findings confirm the equilibrium ensemble equivalence 
01 solely on the basis of Hamiltonian dynamics. For 
close to equilibrium initial conditions, the Nose-Hoover 
dynamics reproduces the caloric curve (Fig. as 
well. 

We now turn to the non-equilibrium properties of the 
model by setting for the HMF system out-of-equilibrium 
initial conditions. In particular, at t — t*, we consider a 
delta distribution for the angles {pHMF{d,t*) — S{0) so 
that fn^HMF^^*) ~ 1) ^'^'i ^ uniform distribution for the 
velocities, pHMF{l,t*) = 1/21, I G [-1,1], with [~ 2.03 
i^HMFit*) — 0.69). If e* = 0, we verified the known re- 
sult that for such initial conditions the systern, after a fast 
process, is dynamically trapped into a QSS Ifll Ull 
(Fig. |2|i). The initial {t < t* + l) violent relaxation corre- 
sponds to a quick mixing of the spins in the single-particle 
/j,-space . The QSS is then characterized by rnj^j^p ~ 
{Thmf = 0.38) for M ^ 00 (zero force) and a hfetime 
tx that increases as a power of M |^ . With respect 
to such QSSs, a crucial issue is to see whether they sur- 



vive when the coupling between the HMF and the TB 
is switched on To address this point, we consider 

e* 7^ but keeping the non-equilibrium initial conditions 
described above. The TB-temperature is first fixed at 
To = 0.38. The time dependence of rn^HMF (^''^?>- Et;C) 
suggests that the QSSs indeed exist even in the canoni- 
cal setup and independently of e* (if e* is small enough). 
Denoting by the life time of the QSSs we found that 
tx ~ Af, with rj that tends to zero as e* increases and 
to the microcanonical estimate given in as e* ^ 
(Fig. |31i). Preliminary evidences 0, suggest that tx 
is also influenced by the "surface effect" parameter 7. 
We remark that during the QSS the long-range system 
does not thermalize with the TB. For example, a consis- 
tent change (10%) of Tq does not alter Thmf and even 
the subset of TB-spins in direct contact with the HMF 
model remains at Tq |l3j |. However, energy fluctuations 
are signiflcantly larger than those due to the algorithm 
pre cision {AEhmf/Ehmf ^ i x 10"^ for M = 10^) 
|13|. This distinguishes the canonical QSSs from the mi- 
crocanonical ones. During these QSSs, the HMF model 
are in a partial equilibrium state at a temperature (spe- 
cific kinetic energy) which is not the one of the TB [l3l |. 
Perhaps this is one reason why a Nose-Hoover dynamics 
with the same out-from-equilibrium initial conditions is 
not capable to reproduce the relaxation to equilibrium. 
In fact, we verified ^3 that in such a case a Nose- Hoover 
dynamics displays very strong fluctuations of the dynam- 
ical variables (e.g., Ehmf) that do not decay with time. 
Another important remark is that classical assumptions 
in mesoscopic stochastic equations seem to rule out the 
existence of such canonical QSS. In fact, a stability anal- 
ysis applied to a Fokker-Planck description of the HMF 
model in both ensembles (canonical and microcanonical) 
shows that anomalous velocity PDFs are (neutrally) sta- 
ble only in the microcanonical ensemble [ij . 

The occurrence of canonical QSSs points towards 
an extension of the ensemble equivalence to some as- 
pects of the non equilibrium properties. We find on 
the other hand that there is a substantial microcanon- 
ical/canonical inequivalence in the relaxation to equi- 
librium process that follows the QSS. For example, 
the final equilibrium specific magnetization changes by 
a factor 4 going from the microcanonical (Fig. |3^) 
to the canonical (Figs. EbjC) simulations, indepen- 
dently of e*. A further indication of this inequiva- 
lence is given by the time evolution of the i/-function, 
H{t) = - dl J^"" de p{l, e, t) \vl{p{1, e, t)), in the two 
ensembles. Indeed, if during the microcanonical dynam- 
ics H{t) on the average increases, reaching its maximum 
at equilibrium, in the e 7^ canonical ones it displays a 
maximum during the QSS and then decreases towards the 
To = 0.38 equihbrium value (See Fig. These uncon- 
ventional behaviors with respect to the iJ-theorem are 
consequences of the different relaxation dynamics in the 
two ensembles, i.e. fixed-energy and fixed-temperature 
(full circles in Fig. [2t). 

In summary, we introduced a Hamiltonian canonical 
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FIG. 4: Time dependence of the _ff-function in the micro- 
canonical (e* = 0) and canonical (e* 7^ 0) ensemble for 
M = 10"* and To — 0.38. The dashed line indicate the com- 
mon maximum of H for e* 7^ reached during the QSS. 

setup for long-range interacting systems through a cou- 
pling with a short-range interacting TB. The coupling 
is described by a parameter that can be tuned continu- 
ously to provide a unified description of microcanonical 
and canonical ensembles. By applying our scheme to 
the HMF model we verified its capability to reproduce 
the equilibrium BG statistics in both ensemble (ensem- 
ble equilibrium equivalence). The major feature of this 
setup is that it is based solely on microscopic Hamilto- 
nian dynamics. This is a novelty with respect to previous 
approaches, that allows an unbiased dynamical descrip- 
tion of the non equilibrium properties of such system. 



As a result we found that, if the coupling with the TB is 
weak enough and we start from out-of-equilibrium initial 
conditions, the dynamics reveals the existence of quasi 
stationary states in the canonical ensemble. These QSSs 
are reminiscent of the microcanonical ones ^ ^, 21e J3 
in the sense that, for example, their lifetime diverges with 
the system size M in a power law fashion. On the other 
hand in presence of the TB, the life-time of the QSSs is 
influenced by the parameters controlling the interaction 
between long-range system and TB. This could be use- 
ful for an experimentalist who is willing to enhance or 
hinder the quasi-stationary behavior Q and could also 
be of some importance in the understanding of the dy- 
namical evolution of quasi-stationary structures, e.g., in 
galaxies 0, or in other long-range interacting systems. 
The presence of canonical QSS extends the notion of en- 
semble equivalence from equilibrium to some non equilib- 
rium properties. A substantial microcanonical/canonical 
inequivalence is found in the relaxation to equilibrium 
process following the QSS and it is clearly revealed by 
a dramatic change in the time dependence of the Boltz- 
mann _ff-function. Of course, a more detailed statistical 
description (and interpretation) of the canonical QSS is 
needed and we are confident that our unbiased set up will 
be a useful tool for this achievement not only with re- 
spect to the HMF model, but also to other Hamiltonian 
long-ran ge s ystems exhibiting either dynamical peculiar 
features |16| or equilibrium ensemble inequivalence Q . 
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